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Abstract 

Using elementary duality properties of positive semidefinite moment matrices 
and polynomial sum-of-squares decompositions, we prove that the convex hull of 
rationally parameterized algebraic varieties is semidefinite representable (that is, 
it can be represented as a projection of an affine section of the cone of positive 
■ semidefinite matrices) in the case of (a) curves; (b) hypersurfaces parameterized 

by quadratics; and (c) hypersurfaces parameterized by bivariate quartics; all in an 
fN| , ambient space of arbitrary dimension. 
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1 Introduction 

ON ' 



Semidefinite programming, a versatile extension of linear programming to the convex cone 
of positive semidefinite matrices (semidefinite cone for short), has found many applica- 
tions in various areas of applied mathematics and engineering, especially in combinatorial 
optimization, structural mechanics and systems control. For example, semidefinite pro- 
gramming was used in [5] to derive linear matrix inequality (LMI) convex inner approxi- 
mations of the non-convex semi-algebraic stability region, and in [6] to derive a hierarchy 
of embedded convex LMI outer approximations of non-convex semi-algebraic sets arising 
in control problems. 

It is easy to prove that affine sections and projections of the semidefinite cone are convex 
semi-algebraic sets, but it is still unknown whether all convex semi-algebraic sets can be 
modeled like this, or in other words, whether all convex semi-algebraic sets are semidefinite 
representable. Following the development of polynomial-time interior-point algorithms to 
solve semidefinite programs, a long list of semidefinite representable semi-algebraic sets 
and convex hulls was initiated in [9] and completed in pp. Latest achievements in the 
field are reported in [7] and [3]. 
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In this paper we aim at enlarging the class of semi-algebraic sets whose convex hulls 
are explicitly semidefinite representable. Using elementary duality properties of positive 
semidefinite moment matrices and polynomial sum-of-squares decompositions - nicely re- 
cently surveyed in [8] - we prove that the convex hull of rationally parameterized algebraic 
varieties is explicitly semidefinite representable in the case of (a) curves; (b) hypersurfaces 
parameterized by quadratics; and (c) hypersurfaces parameterized by bivariate quartics; 
all in an ambient space of arbitrary dimension. 

Rationally parameterized surfaces arise often in engineering, and especially in computer- 
aided design (CAD). For example, the CATIA (Computer Aided Three-dimensional In- 
teractive Application) software, developed since 1981 by the French company Dassault 
Systemes, uses rationally parameterized surfaces as its core 3D surface representation. 
CATIA was originally used to develop Dassault's Mirage fighter jet for the French air- 
force, and then it was adopted in aerospace, automotive, shipbuilding, and other indus- 
tries. For example, Airbus aircrafts are designed in Toulouse with the help of CATIA, 
and architect Frank Gehry has used the software to design his curvilinear buildings, like 
the Guggenheim Museum in Bilbao or the Dancing House in Prague, near the Charles 
Square buildings of the Czech Technical University. 



2 Notations and definitions 

Let P m denote the projective real plane of dimension m, where each element x — [xq, x\, ■ ■ ■ , x m ] ^ 
belongs to an equivalence class [xi/xq, ■ ■ ■ , x m /xo] in M m , the afline real plane of dimen- 
sion m, with Xq = representing the hyperplane at infinity. Let P™ denote the oriented 
projective real plane with non-negative first entry Xq. See [13] for an introduction to 
oriented projective geometry, and its connection with standard projective geometry. 

Let M[x] denote the ring of forms (homogeneous polynomials) in variables x £ P m , with 
coefficients in K, and 

Cd{%) = [^0' ^0 ^1) X ' ' ' ' X X li ' ' " > X m] ^ P ^ ' [%] 

denote a basis vector of m-variate forms of degree d, with s(m, d) = (m + d)\/(m\d\) — 1. 
Let y = [y a ]\a\<2d P s(m > M ) be a real- valued sequence indexed in basis (2d(x), with a £ N m 
and \a\ = a^. A form x i— > p(x) = p T (,2d{%) is expressed in this basis via its coefficient 
vector p £ p s ( m > 2c 0. Given a sequence y £ p s ( m ' 2d ) ; define the linear mapping p i— > L y (j>) = 
p T y, and the moment matrix Mj(y) satisfying the relation L y (pq) = p T Md{y)q for all 
£>,<?£ P s(m ' d) - It has entries [M d (y)]*,P = L y ([Q(x)Ux) T ] a ,is) = y a+p for all at, £ N m , 
+ |/3 1 < 2d. For example, when m = 2 and d = 2 (trivariate quartics) we have 
s(m, 2d) = 14. To the form p(x) = Xq — xqX\x\ + 5x\ X2 we associate the linear mapping 
L y (p) = Uoo — 2/12 + %3i- The 6-by-6 moment matrix is given by 



" 2/oo 
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* 
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2/io 


2/20 
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2/01 


2/n 


2/02 


* 
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* 


2/20 


2/30 


2/21 


2/40 


* 


* 


2/n 


2/21 


2/12 


2/31 


2/22 


* 


. 2/02 


2/12 


2/03 


2/22 


2/13 


2/04 . 
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where symmetric entries are denoted by stars. See [H] for more details on these notations 
and constructions. 

Finally, given a set Z, let convZ denote its convex hull, the smallest convex set containing 
Z. 

3 Convex hulls and moment matrices 

Let 

Z rn4 = conv {C 2d (x) E Pf 1 ' 2 ^ : x E P™} 

and 

y m , d = {ye F s } m > 2d) : M d (y) y o}. 

Theorem 1 Ifm=lord=lord = m = 2 then Z md = y m . d - 

Proof: The inclusion Z m ^ d C y m ,d follows from the definition of a moment matrix since 

M d (( 2d (x)) = Cd(x)Ux f y o. 

The converse inclusion is shown by contradiction. Assume that y* £ Z m d and hence 
that there exists a (strictly separating) hyperplane {y : p(y) = 0} such that p T y* < 
and p T y > for all y E Z m ^ d . It follows that polynomial x h-> p{x) = p T (2d(x) is 
globally non-negative. Since m = 1 or d = 1 or d = m = 2, the polynomial can 
be expressed as a sum of squares of polynomials [H Theorem 3.4] and we can write 
P( x ) = J2 k 1k( x ) = Efc^fcCiO)) 2 = (d(x) T P(d(x) for some matrix P = J^klkQk h 0. 
Then L y {p) = p T y = trace (PM d (y)) = J2k<lkMd(y)qk- Since L y *(p) < 0, there must 
be an index k such that q r [M d (y*)qk < and hence matrix M d (y*) cannot be positive 
semidefinite, which proves that y* ^ y m ,d- d 

In the above definitions we used projective spaces, which have the property of being 
compact, hence closed. If we use affine spaces, Theorem [1] is not correct in this form. 
Indeed, when m = 1 and d = 2, the sequence y* = [1, 1, 1, 1, 2] belongs to y m>d but not to 
Z m ,di as recalled in [SJ Example 5.10]. Theorem [1] becomes correct however by replacing 
the convex hull by its closure in the definition of set Z m ^ d . 

4 Rational varieties 

In projective space P n , a rational variety is the image of P m through a polynomial mapping. 
Define the linear map A : y E p s ( m > M ) i— >• Ay E P n characterized by a matrix A E 
]gmxs(ra,2d)^ ra ^j ona i variety is defined as 

V m4 = A {( 2d (x) E F (m ' M) : x E P m } = {A( 2d (x) E F n : x E P m }. 

Theorem [1] identifies the cases when the convex hull of this rational variety is exactly 
semidefinite represent able. That is, when it can be formulated as the projection of a 
linear section of the semidefinite cone. 
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Corollary 1 If m = 1 or d = 1 or d = m = 2 then 

conv V m , d = {Ay G P» : M d (y) h 0, y G pf 1 ^}. 

Proof: Since conv V m ,d = conv [A((2d(x))) = *A(conv (C2d(^))) = A(Z m> d), the result 
readily follows from Theorem [TJ □ 

The case m = 1 corresponds to rational curves. The case d = 1 corresponds to quadrati- 
cally parameterized rational hyper surf aces. The case d = m = 2 corresponds to hypersur- 
faces parameterized by bivariate quartics. All these rational varieties live in an ambient 
space of arbitrary dimension n > m. 

In all other cases, the inclusion conv V m ,d C A[Z m ^ is strict. For example, when d = 
3, m — 2, the vector G P+ 7 with non-zero entries 

2/00 = 32, |/20 = 2/02 = 34 > 2/40 = 2/04 = 43 > 2/22 = 30 > 2/60 = 2/06 = 128 > 2/4*2 = 2/24 = 28 

is such that M 3 (y*) y but L*(p*) < for the Motzkin form 

x\x\ which is globally non-negative. In other words, y* G A(Z m d ) but y* conv V m ,d- 

5 Examples 

5.1 Trefoil knot 

The space trigonometric curve 

V = {v G IR 3 : ui(a) = cos a+2 cos 2a, V2(a) = sin a+2 sin 2a, v 3 (a) = 2 sin 3a, a G [0, 2tc}} 

is called a trefoil knot, see [2] and Figured] 
Using the standard change of variables 

2x$x\ 

cos a = — ^ 7T, sin a = — t, - 

and trigonometric formulas, the space curve admits a rational representation 

V = {v G P 3 : v (x) = (x 2 + xj) 3 , Vi(x) = (xl + xl)(3x^ - 12x 2 xl + x\), 

v 2 {x) = 2x a;i(xo + xf)(5xl — 3xf), v 3 (x) = Axox^xl — 3xf)(3xl — xf), x G P} 

as the image of the projective line P through a sextic mapping, i.e. n = 3, m = 1 and 
d = 3 in the notations of the previous section. 

By Corollary [T], the convex hull of the trefoil knot curve is exactly semidefinite repre- 
sent able as 

conv V = {Ay G P 3 : M 3 (y) >z 0, y G P+} 

with 



A = 



1 





3 





3 
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-11 
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12 
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Figure 1: Tube plot of the trefoil knot curve, whose convex hull is exactly semidefinite 
representable with 3 liftings. 



and 



M 3 (y) 



yo * 

yi 2/2 

2/2 2/3 

2/3 2/4 



* 
* 

in 
ys 



* 
* 

ye 



where symmetric entries are denoted by stars. The linear system of equations v 
be solved by Gaussian elimination to yield the equivalent affine formulation: 



Ay can 



conv V = {v E M 3 : 

|(3 + vi + 2 Ul - Au 3 ) 
^{10v 2 + v 3 + 48m 2 ) 
^(3-v 1 -20w 1 -2« 3 ) 



2^(6w 2 - 5w 3 + 96w 2 ) 



Ui 



2^(6f 2 - 5v 3 + 96u 2 
which is an explicit semidefinite representation with 3 liftings 



* 
* 

"2 



* 
* 
* 

U 3 



y o, u e M 3 } 



5.2 Steiner's Roman surface 

Quadratically parameterizable rational surfaces are classified in [3]. A well-known example 
is Steiner's Roman surface, a non-orientable quartic surface with three double lines, which 
is parameterized as follows: 

V = {v e R 3 : v 1 {x) = 2 , v 2 (x) = 2 , v 3 (x) = .^ X Tl 2 , x e M 2 } 

1 + x{ + x\ 1 + x{ + x\ 1 + x\ + x 2 

see Figure [21 



Figure 2: Two views of Steiner's Roman surface, whose convex hull is semidefinite repre- 
sentable with 2 liftings. 



In projective coordinates, the surface becomes 



V = {v G P 3 : Vo(x) = xl + xl+xl, V\(x) = 2x Xi, v 2 (x) = 2x x 2 , v$(x) = 2x x x 2 , x G P 2 } 

which the image of the projective plane P 2 through a quadratic mapping, i.e. n = 3, 
m = 2 and d = 1 in the notations of the previous section. By Corollary [U its convex hull 
is exactly semidefinite representable as 



conv V = {Ay G P 3 : M^y) h 0, y G P^} 



with 



^4 



and 



MM 



10 10 1 

2 

2 

2 



2/oo * * 
Z/io 2/20 * 
2/oi 2/n 2/02 



The linear system of equations v = Ay can easily be solved to yield the equivalent affine 
formulation: 



conv V = {v G M 3 
which is an explicit semidefinite representation with 2 liftings. 



1 — Mi — u 2 



y 0, u G M 2 } 



5.3 Cay ley cubic surface 



Steiner's Roman surface, studied in the previous paragraph, is projectively dual to Cay- 
ley's cubic surface {v G P 3 : detC(v) = 0} where 



C(v) 



Vq * * 
V\ Vq * 
V 2 V 3 V 



The origin belongs to a set delimited by a convex connected component of this surface, 
admitting the following affine trigonometric parameterization: 

{v G 1R 3 : vi(a) = cosai, v 2 (a) = sina 2 , 

1*3(0;) = cosai sina 2 — coso^sinai, aci G [0, tt], a 2 G [— 7r,7r]}. 



V 



This is the boundary of the LMI region 

conv V = {veF 3 + : C(v) y 0} 

which is therefore semidefinite representable with no liftings. This set is sometimes called 
a spectrahedron, a smoothened tetrahedron with four singular points, see Figure |3j 



Figure 3: Convex connected component of Cayley's cubic surface, semidefinite repre- 
sentable with no liftings. 



Using the standard change of variables 

Xq X; . 2x$Xi . 

cosaj = — ^ Tt, sin«j = — ^ k, % = 1,2 

we obtain an equivalent rational parameterization 

V = {v G P 3 : Wo (a;) = (^o + ^i)( x o + x f)i v i( x ) = ( x o ~ x l)( x o + x l)-> 

v 2 (x) = 2x x 2 (xl + x\), v 3 (x) = 2x (-xi + x 2 )(xq + xix 2 ), x G P 2 }. 
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which is the image of the projective plane P 2 through a quartic mapping, i.e. n — 3, 
m = 2 and d = 2 in the notations of the previous section. By Corollary (TJ its convex hull 
is exactly semidefinite representable as 

conv V = {Ay eF 3 + : M 2 (y) y 0, y e P+ } 

with A of size 4-by-15 and M 2 (y) of size 6-by-6, not displayed here. It follows that conv V 
is semidefinite representable as a 6-by-6 LMI with 11 liftings. 

We have seen however that conv V is also semidefinite representable as a 3-by-3 LMI with 
no liftings, a considerable simplification. It would be interesting to design an algorithm 
simplifying a given semidefinite representation, lowering the size of the matrix and the 
number of variables. As far as we know, no such algorithm exists at this date. 

6 Conclusion 

The well-known equivalence between polynomial non-negativity and existence of a sum- 
of-squares decomposition was used, jointly with semidefinite programming duality, to 
identify the cases for which the convex hull of a rationally parameterized variety is exactly 
semidefinite representable. Practically speaking, this means that optimization of a linear 
function over such varieties is equivalent to semidefinite programming, at the price of 
introducing a certain number of lifting variables. 

If the problem of detecting whether a plane algebraic curve is rationally parameterizable, 
and finding explicitly such a parametrization, is reasonably well understood from the 
theoretical and numerical point of view - see [H] and M. Van Hoeij's algcurves Maple 
package for an implementation - the case of surfaces is much more difficult [12]. Up to 
our knowledge, there is currently no working computer implementation of a parametriza- 
tion algorithm for surfaces. Since an explicit parametrization is required for an explicit 
semidefinite representation of the convex hull of varieties, the general case of algebraic 
varieties given in implicit form (i.e. as a polynomial equation), remains largely open. 

Finally, we expect that these semidefinite representability results may have applications 
when studying non-convex semi-algebraic sets and varieties arising from stability condi- 
tions in systems control, in the spirit of [3 [6] . These developments are however out of 
the scope of the present paper. 
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